def tau(y):    return alpha*(1-y)/2 + betta*(1+y)/2
def ro(y):    return ro1(tau(y))
def z(y):    return z1(tau(y))
def dro(y):    return dro1(tau(y))*mult
def dz(y):    return dz1(tau(y))*mult
def ddro(y):    return ddro(tau(y))*(mult**2)
def ddz(y):    return ddz1(tau(y))*(mult**2)
def dddro(y):    return dddro1(tau(y))*(mult**3)
def dddz(y):    return dddz1(tau(y))*(mult**3)



m = 10.
l = arange(2*m+1)
phi = 2*pi/(2*m+1)*l
phi[0] = 0.0001
print f_Smooth_Im(2,5,phi).sum()*2*pi/(2*m+1)

print K1_tau(1,1,0)
xx, yy = meshgrid(ind1, ind1)
xx[diag_indices_from(xx)] = 0
yy[diag_indices_from(yy)] = 3
from time import time
t1 = time()
KK = K1_tau(1,yy,xx)
print KK
t2 = time()
print t2-t1